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ABSTRACT 

We present evidence for cosmofogical gas accretion onto spiral galaxies in the local 
universe. The accretion is seen through its effects on the dynamics of the extra-planar 
neutral gas. The accretion rates that we estimate for two nearby spiral galaxies are of 
the order of their star formation rates. Our model shows that most of the extra-planar 
gas is produced by supernova feedback (galactic fountain) and only 10 — 20% comes 
from accretion. The accreting material must have low specific angular momentum 
about the disc's spin axis, although the magnitude of the specific angular-momentum 
vector can be higher. We also explore the effects of a hot corona on the dynamics of 
the extra-planar gas and find that it is unlikely to be responsible for the observed 
kinematical pattern and the source of accreted gas. However, the interaction with the 
fountain flow should profoundly affect the hydrodynamics of the corona. 

Key words: galaxies: kinematics and dynamics galaxies: individual: NGC 891, NGC 
2403 galaxies: haloes galaxies: evolution ISM: kinematics and dynamics 



1 INTRODUCTION 



several spiral galaxies have revealed 
or HI haloes surrounding the galac- 



'Swaters, Sancisi & van dor Hulstl 1 19971 : 



Observations of 
> thick HI layers 
' tic discs (e.g . 

' iFraternali et al.l [2b02: Matthews fc Wood .200a '). The kine 
matics of this extra-planar gas is characterized by: 
■ (i) a decrease in rotation velocity in th e ver- 

' tical direction llSwaters. Sancisi fc van der Hulstl 1 19971 : 
IFraternali et al.l 12005); (ii) vertical motions from and 



towards the disc jBoomsrnal 20071 ): (iii) a general 
dial inflow (IFraternali et al.l l200ll ). In a previous paper 



l|Fraternali &: Binnev 20061 . hereafter FB06) we developed 
the galactic fountain model of extra-planar H I in which gas 
clouds are shot upwards from regions of star formation in the 
disc, and then follo w ballistic trajectories before re turning 
to the disc (see also lCollins. Benjamin fc Randll2002l ). Com- 
paring the predictions of this model with data from deep H I 
surveys of the edge-on galaxy NGC 891 and the moderately 
inclined galaxy NGC 2403 we concluded that the model fails 
in two respects: (i) it under-predicts the difference between 
the rotation rates of disc and extra-planar gas, and (ii) it pre- 
dicts that the extra-planar gas should on average be flowing 
out rather than in. We showed that these failures could not 
be rectified by assuming that the clouds were visible in the 
21-cm line on only part of their trajectories. 

In this paper we extend the fountain model to include 
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interaction with extragalactic gas. This gas is expected to 
take two forms: (i) the hot corona of gas at the virial tem- 
perature, and (ii) cold, infalling clouds. Standard cosmolog- 
ical models of galaxy formation imply the existence of both 
types of gas. In fact, arguments based on the cosmic back- 
ground radiation and primordial nucleosynthesis, combined 
with the measured baryonic masses of galaxies imply that 
most of the baryons in the Universe must lie in ext ragalactic 
gas (e.g. lWhite fc Frenkll99ll : ISommer-Larsen|[2009' ). Limits 
on the background radiation levels at frequencies from the 
X-ray to radio bands imply that most of the gas must be 
at low densities and hot, T ^ 10^ K. From the existence of 
Ca II absorpt ion lines in the spectra of stars at high Galac- 
tic latitudes, ISpitzerj (|l956l ) inferred the existence of such 
gas around the Milky Way, and by mappin g the soft X-ray 
emission from the massive galaxy NGC 5746 Pedersen et al.l 
(2006) have establis hed the existence of such gas around at 
least some galaxies. iBeniamin fc Danlvl (|l997l ) showed that 
observed H I clouds of fountain gas will experience apprecia- 
ble drag as they move through the coronal gas. The head-tail 
morphology of many co mpact high-velocity clouds around 
the Milky Way (|Bruns et al. 2001) confirms this prediction. 
In this paper we add this drag to our earlier dynamical model 
and show that although it cannot account for the observed 
velocities of extra-planar H I, it must have major implica- 
tions for the hydrodynamics of the coronal gas. 

The existence of streams of cold gas around galaxies is 
less securely established, but is supported by several lines of 
argument. From a theoretical perspective, extended Press- 
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Schechter theory l|Bond et al.lll99ll : lLacev fc Cole|[l993l ') pre- 
dicts that galaxies grow through a series of infall events. 
These range from a small number of major mergers, right 
down to an almost continuous drizzle of infalling dwarf 
galaxies and gas clouds. It is hard to assess the impor- 
tance of this drizzle from numerical simulations of limited 
mass and spatial resolution, but as simulations have become 
more sophisticate d, awareness of the im portance of the driz- 
zle has increased llBirnboim fc Dekej2003.:,Keres et al.„2005i : 



ISemehn fc Combesll2005l : iDekel fc Birnboimll200 

The rate of star-formation in the solar neighbo rhood 
has been remarkably constant over the Galaxy's life (T warogj 
ll980l : lBinnev. Dehnen fc Bertem1l2000l l. which suggests that 
the gas density has not decreased significantly even though 
the current star-formation rate would exhaust the gas in 
a couple of gigayears. These facts strongly suggest that gas 
consumed by star formation is mostly replaced by accretion, 
whether by cooling of the corona or cold infall. Steady ac- 
cretion of metal-poor gas would also explain the discrepancy 
between the observed stellar metallicity distribution in the 
solar neighbourhood and that^predicted by the closed-box 
model of chemical evolution (Tinslev. .l981 ; Mattcucci 2003). 

What remains unclear is whether gas is predominantly 
accreted from the hot corona or comes by cold infall, and 
in the latter case, how important is steady drizzle compared 
to the acquisition of gas in parcels associated with discrete 
events such as the capture of the Saggitarius dwarf galaxy. 
In our Galaxy, a limit can be placed on the masses of in- 
falling o bjects because the solar nei ghbourhood, which is 
both old (|Binnev. Dehnen fc Bertelli I i2000) and cold, would 
be rapidly heated by the passage through it of objects with 
masses > lO'^ Mq (|Lacev fc Ostrike j| 19851 : iToth fc Ostrikeij 
1 19921 ). 

An observational indication of the importance of cold 
infall is the fact that the Galaxy is surrounded by a 
large nu mber of clouds, the so-called h igh-velocity clouds 
(HVCs) IWakker fc van WoerdenI ll997D that can be ob- 
served in H I and on the average have negative radial veloc- 
ities. New distance estimates for the large HVC complexes 
show that they are located in the upper hal o of the Milky 
Way and ha ve masses of order lO'^^'^ Mq l|Wakker et all 
l2007al . 120081 '). Most HVCs are too far from the plane to 
be part of the Galactic fountain, and those that have 
measured metallicities are too metal-poor to be predomi- 
nantl y gas ejected from the star-forming disc (|Tripp et al.l 
I2OO3I ). The net infall rate M associat ed with high-velocity 
clouds seems to be M ~ 0.2 Mq vr~^ JWakker et al.|[2007al : 



I 

iPeek. Putman fc Sommer-LarsenI 120071 ). which is about an 

order of magnitude lower than the SFR of the Milky Way. 

There are observations of infalling g as complexes 
also in external galaxies (for a review see ISancisi et al.l 
|2008| ). A well known example is the m assive HI com- 
plex accreted by M 101 (|van der Hulst fc: Sancisi.,19 88) but 
the majority of these accreting clouds will have much 
lower masses and have escaped detection. In fact, new 
deep observations of the region of space around M31 
reveals the presence of nu merous HVCs with masses 
down to a few 10'' Mq (|Westmeier. Braun. fc Thilkerl 
I2OO5I ). In some cases, halo clouds are detected at 
very anomalous [counter-rotating) velocities with re- 
spect to the disc rotation and massive cold filaments 
are also observed l|Oosterloo. Fraternali fc Sancis 1 I2OO7I : 



iFraternali et al.ll2002l '). Gas accretio n may be also arise from 
interactions with satellit e galaxies l|van der Hulst fc Sancisil 
I2OO5I : ISancisi et al1l2008l ') . 

The above data suggest that star forming galaxies are 
accreting material from th e intergalactic spac e at a rate of 
at least M ~ 0.2 Mq yr"^ (|Sancisi et al.ll2008l ') which is typ- 
ically 10 — 20% of the SFR for these objects. In this paper 
we argue that this directly observed accretion rate is in fact 
a lower limit. Most of the accreting gas interacts with gas 
pushed up into the halo by stellar feedback, and it is 06- 
servahle only indirectly via the peculiar kinematics of the 
extra-planar gas. 



2 THE MODEL 



Our p ure fountain model is described in lFraternali fc BinnevI 
(|2006l '). We integrate orbits of particles that are ejected from 
the disc of a spiral galaxy and travel through the halo until 
they fall back to the disc. The data require that the particles 
are ejected nearly straight upwards. The kinetic energy of 
ejection is ^ 4 percent of the mechanical energy produced 
by supernovae. The potential of the galaxy comprises the 
contributions of four components: stellar and gaseous discs, 
both described by exponential profiles, a dark matter (DM) 
halo and a bulge, both descri bed by double power law pro- 
files (|Dehnen fc Binnevlll998|). The potential is constrained 
by the HI rotation curve jOosterloo. Fraternali fc Sancisil 
I2OO7I : iFraternaU et al]|2002l '). After each timestep of the or- 
bit integrations, the positions and velocities of the particles 
are projected along the line of sight to produce an artifi- 
cial cube of positions on the sky and radial velocities. This 
cube is then smoothed and compared with the observed H I 
data cubes. Typically, the mass and spatial resolutions of the 
simulations are two orders of magnitude higher than those 
of the data, therefore the global effect is that of a smooth 
fountain flow (see also FB06). 



3 INTERACTION WITH THE HOT CORONA 

In this section we concentrate on the edge-on galaxy 
NGC 891, for which parameters of the corona can be ob- 
tained from X-ray ob servations l|Bregman fc Houckl [l997l : 
IStrickland et al.ll2004l ). We model the corona as a plasma 
of temperature T that is in equilibrium in the axisymmetric 
potential of the galaxy $(-R, z). Then the hydrostatic equa- 
tions in cylindrical coordinates imply that the density is 



A exp 



/imp 
' kT 



(1) 



where A is a constant, /i = 0.61 is the molecular weight 
of the plasma, and vo is the speed at which the corona ro- 
tates (assumed to be independent of radius). We adopt the 
tempe rature T = 2.7 x 10^ K determined from X-ray spec- 
tra by [Strickland et al.l (|2004| ) and for several values of vo 
determine A by fitting the predicted X-ray luminosity Lx 
to that observed. Lx was obtained by integrating the emis- 
sion predicted by equation ((TJ over roughly the same re- 
gions as that of the X-ray obse rvations. The emissivity was 
taken from the MEKAL recipe (jMewe. Lemen fc Oordlll986l : 
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iLiedahl. Osterheld k, Goldstein|[l995h that is embedded in 
the XSPEC package. 

The X-ray distribution is not highly flattened, so the 
adopted values of vq lie significantly below the local circular 
speed. 

In reality a significant fraction of the observed X-ray lu- 
minosity will arise from interfaces between cold clouds and 
the corona rather from the main body of the corona. Conse- 
quently, our procedure for determining A will tend to over- 
estimate the corona's density. This fact will reinforce the 
conclusion reached below, that the corona cannot absorb 
angular momentum fast enough to account for the observed 
rotational lag of the H I halo. 



3.1 Cloud dynamics 

We consider the motion of an HI cloud with temperature 
< 10"* K through the corona. Near the disc, the corona's 
cooling time fcooi ~ mpkT/{Ap) is }^ 100 Myr, which is 
much longer than the flow time 



D/v ~ IMyr 



(2) 



around a cloud of diameter D ~ 100 pc at speed v ~ 
lOOkms"^, so most of the coronal gas is expected to pass 
right by, and we should idealize the interaction as the flow 
of a gas around a dense object. The densities of the cloud 
and the corona differ by a factor ^ 300, so the motion is 
akin to the flight through air of water drops from a garden 
sprinkler. 

We assume that the flow of the corona past a cloud has a 
high Reynolds number; this assumption minimizes the drag 
on the cloud, and we shall see that even this minimum drag 
is highly significant. When the Reynolds number is high, the 
drag is proportional to the square of the flow velocity v of 
the corona past the cloud, so the acceleration of the cloud 
can be written 



(3) 



where C is a dimensionless constant of order unity, and L 
is the distance over which the cloud interacts with its own 
mass of coronal gas: 



m 
pa 



(4) 



with a and m the cloud's cross-sectional area and mass, 
respectively. In the absence of other forces, drag causes the 
cloud's velocity to decay as 



Vo 



1 + Ctvo/L' 



(5) 



Thus the speed halves in a time tdrag = L/{Cvo) ~ SOOtflow. 
For t;o ~ 100 km s"'^, tdrag ~ 300 Myr is of the same order 
as the orbital times of the clouds that make up the H I halo 
(FB06 Fig. 10). For low Reynolds numbers, tdrag would be 
even lower than the above estimate. Note that if a cloud 
were to sweep up all the coronal gas that it encounters rather 
than deflecting it, the drag on a cloud would be still given 
by equations ((3} and Q. 




Figure 1. Azimuthal velocities acquired by the hot corona per 
Myr due to the transferring of angular momentum from the foun- 
tain gas. From outside in contours are at 0.25, 0.5, 1, 2, 4, 8, 16, 
32, 64 kms~VMyr. 



3.2 Coronal spinup 

Since these numbers establish that momentum transfer be- 
tween clouds and coronal gas may be a significant process, 
consider the possibility that this process produces the addi- 
tional lag of the H I halo relative to the disc that FB06 found 
the observations to require. Let fountain gas be ejected at a 
rate Mf and let S£ be the amount by which its specific an- 
gular momentum is reduced by interaction with the corona. 
Then the rate at which the corona absorbs angular momen- 
tum is L = M{5i, and the time required for the corona to 
be spun up to near the circular speed is 



A:/corona^c 



(6) 



where Ic — Rvc is the specific angular momentum of a cir- 
cular orbit. 

Considering just the HI that is more than 1.3 kpc 
from the plane, FB06 found that S£/£c ~ 0.05 and Mi = 
lO^M0Myr-\ so t.pin ~ (2Afcorona/10^ M©) Myr. The 
mass of the corona inside radius r is 3 x lO^(r/lOkpc)M0, so 
tspin ~ 60(r-/10 kpc) Myr. That is, the capacity of the corona 
for absorbing angular momentum is so slight that it will be 
brought to corotation with the HI halo within a dynamical 
time unless the angular momentum is moved from the radii 
~ 10 kpc at which it is imparted to much larger radii. This 
cannot be done in less than a dynamical time. 

To test this conclusion, we added the deceleration ^ 
to the fountain model presented in FB06 - a description of 
the initial parameters of the model can be found in FB06. 
Here the potential is taken to be the one appropriate to a 
maximum disc, but in FB06 we showed that the kinematics 
of the extra-planar gas is insensitive to the choice of poten- 
tial. As determined in FB06, the total mass of the H I halo 
is Afhaio = 2 X lO'^M©. 

We adjusted the ratio m/a until the model reproduced 
the observed lag velocity. We considered hot coronas rotat- 
ing at difi^erent speeds and with different velocity patterns. 
In the case of a corona that rotates at lOOkms"^ indepen- 
dent of radius, the data for NGC 891 are s uccessfully simu- 
lated when m/a ~ 10^ Mq kpc~^ (see also iFraternali et al.l 
l2007f) . 

However, in this picture, in which the coronal gas simply 
flows around clouds without being accreted onto them, all 
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the angular momentum lost by the H I clouds of the halo has 
to be absorbed by the corona. If angular momentum taken 
up by the corona in an annulus centred on {R, z) remains 
in that annulus, then in 1 Myr the azimuthal velocity of 
the corona at [R, z) will increase by the amounts plotted 
in Fig. [T] we see that the spin-up time of the inner corona 
is extremely short; within a time tspin ~ 10 Myr the inner 
corona will co-rotate with the disc and the drag will vanish. 

Thus the corona can be responsible for the observed 
rotational lag of the HI halo only if the clouds sweep up 
the coronal gas along their path rather than deflecting it, 
which is not expected to happen given the large ratio of the 
cooling time of coronal gas to the flow time around a cloud. 
If clouds do not sweep up coronal gas, and the ratio m/cr is 
not much larger than the above estimate, the HI halo will 
stir the corona rather like the impeller of a centrifugal pump 
stirs the fluid in the pump. In response to this stirring, the 
coronal gas will move rapidly outwards parallel to the disc 
plane, and be replaced by gas that moves down parallel to 
the symmetry axis. Thus the H I halo must drive a large-scale 
circulation in the corona. Calculation of the dynamics of this 
circulation lies beyond the scope of this paper, but is a topic 
that needs to be pursued because it must have significant 
implications for the corona and X-ray observations. 



4 SWEEPING UP OF AMBIENT GAS 

In light of the conclusion above that the hot corona can be 
responsible for the observed rotational lag only if it is swept 
up by fountain clouds, we now examine the case in which 
fountain clouds sweep up ambient gas. Given the long cool- 
ing time of the coronal gas, we envisage that the bulk of 
the gas swept up is cold, infalling material. When a foun- 
tain cloud impacts on a stream of cold gas, some material 
is likely to be heated to the virial temperature, but given 
the high densities involved, and the potential for rapid tur- 
bulent mixing of the two bodies of gas, other material will 
be swept up into the cloud of fountain gas. For simplicity 
we assume that this is the dominant process. That is, we 
assume that a travelling cloud of mass m gains mass at a 
rate m = am and momentum at a rate rhvi, where Vi is 
the velocity of the infalling material. We implement these 
assumptions by modifying the procedure used in FB06 to 
determine the trajectories of clouds as follows. After each 
timestep of duration 5t we increment the cloud's mass by 
5m = amSt and change its velocity from vo to 

movo + Smvi 

'^i = . (7) 

mo -f dm 

where Vi is the velocity of the infalling material. 

The parameter a, which is proportional to the volume 
density of infalling material, becomes the key parameter to 
be constrained by the data. Its spatial variation a(R, z) is 
not known. However, experimenting with different patterns 
we found that the results are not qualitatively very different 
as long as its value does not vary too steeply in the inner 
halo. We therefore take a to be constant, which implies that 
in a given interval 5t a cloud captures the same amount of 
infalling gas regardless of its location. We exclude regions 
very close to the disc, \z\ < 0.2 kpc. 



A more critical characteristic of infalling material is the 
velocity Vi at which it is captured. At any given time infalling 
material is likely to have some non- vanishing average angu- 
lar mo mentum, but the d i rectio n of this average changes over 
time iQuinn fc BinnevI Il992l : lAubert. Pichon fc Colombil 
|2004 ). and we neglect the correlation between its average 
value at the present time and the angular-momentum vec- 
tor of the disc. 

We consider two extreme possibilities: one possibility 
is that the z component of Vi is dominant and negative, 
implying that infall is predominantly parallel to the z axis. 
The modulus of Viz is chosen at random between zero and 
the local escape speed. The other two components of «i are 
taken to follow Gaussian distributions with zero mean and 
dispersion ~ 50kms~^. The magnitude of vi is constrained 
to be less than the escape speed. With this prescription, the 
angular momentum about the z axis of an infalling cloud 
vanishes in the mean. 

The second possibility considered is that the dominant 
component of Vi is parallel to the local radius vector. For this 
component we use the distribution that was previously used 
for Viz . The two tangential components of Vi are drawn from 
a Gaussian distribution with dispersion ~ 50kms~^, again 
subject to the constraint that the total speed should not 
exceed the escape speed. This prescription may be motivated 
by the observation that in cosmological simulat ions infalling 
material is typically on highly eccentric orbits (iKnebe et al.l 
l2004h . On such an orbit, the radial component of the velocity 
vector is usually dominant. We also considered models with 
random velocity components and no preferred directions; 
these gave results in between the two extreme possibilities 
described here. 

4.1 Application to NGC 891 

We applied our model of a fountain that interacts 
with infalling gas to the HI obse rvations of NGC 891 
(jOosterloo. Fraternali fc Sancisi|[2007l ). As in FB06, clouds 
were ejected nearly vertically from the disc with a veloc- 
ity drawn from a Gaussian distribution, and the dispersion 
of this distribution was adjusted to optimize the agreement 
between the model and observational data cubes. 

Fig. [5] (leftmost column) shows five sample channel 
maps for NGC 891 at 28" (about 1.3 kpc) resolution; helio- 
centric velocities in kms~^ are in the top-left corners. A key 
feature of the data is that the apparent vertical thickness of 
the channel maps increases as the velocity moves away from 
values characteristic of fast rotation towards the systemic 
velocity. Thus in the Whci = 292kms~^ channel the HI ap- 
pears thin, and in the Vsya = 528kms~^ channel it has a 
very extended vertical distribution. This behaviour reflects 
the lag of the HI halo with respect to the disc and is the 
main constraint on our models. 

In the rightmost column of the same Fig. [5] we report 
the prediction from a model in which the extra-planar gas is 
produced by stellar feedback only (pure galactic fountain), 
presented in FB06. Note t hat, thanks to the availability 
of a new deeper data cube (|Oosterloo. Fraternali fc Sancisil 
|2003) a lower contour is shown than appears in FB06. These 
upper model maps show the distribution of H I to be much 
thicker than is found in the observed channel maps on the 
extreme left. This mismatch arises because in the model of 
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Figure 2. Comparison between 5 channel maps for NGC 891 jOosterloo. Praternali &: Sancis 1 120071) and those obtained with three 
dynamical models. The first column shows the data, heliocentric radial velocities (in kms~^) are reported in the upper left corners. The 
channel map in the bottom row is roughly at the systemic velocity (Vsys = 528kms~^). Contour levels (for data and models) are: 0.25 
(2.5cr), 5, 1, 2, 5, 10, 20, 50 mjy beam-^ 



Table 1. Fountain + accretion models for NGC 891 



Model 


Phase 
change 


/iv 

(kms-1) 


a 

(Gyr-i) 


Min 


Mout 


Polar 
Radial 


no 

yes 


90 
80 


2.0 
0.6 


3.4 
2.3 


0.12 
0.05 



FB06 the rotation velocity decreases too slowly witli increas- 
ing z. FB06 show that this problem cannot be eliminated by 
using different initial conditions or supposing that the gas is 
invisible in the first part of its orbit (phase-change models) . 

Consider now the models shown in the two central 
columns of Fig. [2] In these models the clouds of the fountain 
interact with an accretion flow as described above, and dif- 
fer from each other in the dominant component of the infall 
velocity Vi: in the left column Vi^ is dominant, while in the 
right column vw dominates. The specific accretion rates are 



Q = 2 Gyr ^ and a = 0.6 Gyr ^ respectively for the two 
models (see Table [1} 

It is evident that the top two channel maps of the left- 
central column have narrower vertical distributions than are 
obtained without infalling gas, and thus provide much better 
fits to the data. Modest improvements in the shape of the 
lower two channel maps for velocities near systemic are also 
evident. We find that the total accretion rate of cold gas 
required by the model with dominant polar component of 
accretion is Minflow — 3.4MQyr~'^, which is very close to 
the star formation rate of NGC891 (SFR = 3.8MQyr"^ 
iPopescu et al]|2004D . 

When the dominant component of the infall velocity is 
radial, the data are best fitted when the fountain gas suffers 
a phase change such that it becomes visible as HI when \z\ 
peaks. A galaxy with the star-formation rate of NGC 891 
produces a significant fiux of ionizing radiation and it is 
plausible that a large fraction of the outgoing fountain fiow 
would be ionised. The right-central column shows channel 
maps for a model in which cold material is accreting at a 
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rate Minflow — 2.3M0yr~ . We see that this model fits the 
data as well as the model with polar accretion shown in the 
left-central column. 

The two accretion models presented here are best fits 
to the data obtained via a minimization of the residuals be- 
tween the data and the models, much as described in FB06 
except that here the discrepancies are also considered chan- 
nel by channel to the whole data cube and not only in the 
total HI maps. Once the mass of the halo has been fixed 
(Afhaio — 2 X 10^ Mq), the key parameters are the accretion 
parameter a (related to the accretion rate Minflow) and, as 
in the pure fountain model, the characteristic kick veloc- 
ity. For the latter we found values slightly higher than in 
the pure fountain model (80— 90 km s~^) because to reach a 
given height fountain gas has to overcome drag from infalling 
gas in addition to the galaxy's gravitational field. The total 
energy input from the disc required to produce such a foun- 
tain fiow is ^6 ( ^ 8 for the phase change model) per cent 
of the energy available from the supernovae. These percent- 
ages are a factor ~ 1.5 larger than those of the pure fountain 
model (FB06). The actual value of the supernova feedback 
in galaxies is uncertain; our valu es agree with the predictio n 
of hydrodynamical models (e.g. iMac Low fc Ferraralll999l ). 
The infiowing (accretion) is only ~ 10 percent of the out- 
flowing (fountain) gas (see Table [l|. 

In the model with radial accretion and phase change 
the accretion rate is not well constrained and the data al- 
low for values of a between 0.2 and 1.5 Gyr~^. This is not 
surprising since this phase-change model produces a fairly 
good fit to the data of NGC891 also without accretion, as 
shown in FB06. However in the case of low accretion rate it 
fails completely to reproduce the inflow pattern observed in 
NGC 2403. We show in the next section that the inclusion 
of accretion elegantly solves this problem. 

In the above models we assumed that the angular mo- 
mentum of the accreting material about the disc's spin axis, 
Lz, is zero on average. We have tested whether this is a nec- 
essary requirement for the model and found that the data 
allow for a non-negligible Lz- At radius R the material ac- 
creting on a disc with circular speed Vc must have on average 
Lz < ^Rvc- If I/z = ^Rvc, the accretion rate required to re- 
produce the data becomes Minflow ~ 6 Mq yr~^ . We discuss 
the implications of this requirement in Section [G] 

iFraternali et al.l (|2005l ) determined the rotation speed 
Vrot of gas in NGC 891 at several heights above the plane 
by identifying the terminal velocity VteTin{R,z) at each po- 
sition {R,z). The rotation speed is taken to be Vrot ~ 
Vterm — v^cr?an + o-Qbs ' whcre CTran IS the raudom velocity of 
HI clouds and o" oba = 6kms~^ is the observational error. 
IFraternali et al. I l|2005l ) assumed that aran — 8kms ev- 
erywhere. More recent modelling shows clear evidence for 
a larger dispersion gran — 20kms~^ away from the plane 
jOosterloo. Fraternali fc San cisill200 7|). So i t is a ppropriate 
to move the points from Fraternali et al.l l|2005h down by 
llkms"^ = VW+W - (details in Fraternah, in 

prep.); the data points in Fig.|3]show the corrected points for 
z = 3.9 and z = 5.2 kpc. The lines show model predictions 
obtained from weighted means of the azimuthal components 
of the particle velocities. We see that the model with infall 
(full/blue curves) fits the data points very well, while the 



250 



200 



150 



o 



B 
o 



100 



50 



'^'^'^ Fountain 
Fountain + Accretion 




z=3.9 l<pc 



2 4 6 8 10 12 14 16 18 
R (kpc) 



250 



200 



150 

o 
o 

03 

c 100 
o 

o 

50 



Disc . 



Fountain 



Fountain + Accretion 




z=5.2 l^pc 



2 4 6 8 10 12 14 16 18 
R (kpc) 

Figure 3. Rotational velocities (points) obtained from the HI 
data at 3.9 kpc (u pper panel) and 5.2 kp c (lower panel) from the 
plane of NGC 891 (praternali et al.ll2005h . The dotted line shows 
the rotation curve in the plane, while the other two lines show 
the azimuthal velocity predicted by our pure fountain (dotted- 
dashed/red) and fountain+accretion (solid/blue) models. 



pure fountain model (dashed/red curves) yields only half 
the lag required by the data. 



4.2 Application to NGC 2403 

Consider now the ability of the model that includes infall 
to reprodu ce the data for the m oderately inclined galaxy 
NGC 2403. iFraternah et al.1 (|2002l ') showed that the H I mea- 
surements of this galaxy reveal an extended layer of extra- 
planar gas that rotates less rapidly than the disc gas and 
fiows in towards the centre of the galaxy. In FB06 we showed 
that pure fountain models in which the halo gas rotates less 
fast than the disc show outfiow rather than inflow. Here 
we show that adding accretion similar to that derived for 
NGC 891 causes the model to predict inflow similar to that 
observed in NGC 2403. 

Fig. |4]shows position- velocity plots along the minor axis 
of NGC 2403 (middle row) and parallel to this axis at offsets 
of ±2'. The leftmost column shows the data, the rightmost 
column the pure fountain model and the middle column the 
fountain with accretion. Similar results are obtained regard- 
less of whether the accretion has a polar or radial pattern, 
and in the latter case a phase change is not required (see Ta- 
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Figure 4. Comparison between 3 p-V plots parallel to the minor 
axis of NGC 2403 and those obtained with two dynamical models. 
The first column shows the data, contour levels (for data and 
models) are 0.5 (~2.5o-), 1, 2, 5, 10, 20, 50 mjy beam"!. 



Table 2. Fountain+accrction models for NGC 2403 



Model 


Phase 
change 


(kms-l) 


a 


Minflow 

(Moyr-i) 


Min 
Mout 


Polar 


no 


55 


2.0 


1.0 


0.21 


Radial 


no 


50 


1.2 


0.6 


0.15 



ble[2]). The results shown are obtained for the same values of 
the accretion parameter a as was derived from NGC 891 (for 
the radial infall model we doubled the value of a since clouds 
are now visible along the whole trajectory). We see that the 
model with infall reproduces the data quite well. The re- 
quired accretion rate is Minflow ~ 0.8 Mq y r~^. For compar- 
ison the SFR in NGC 2403 is 1.2 Mq yr"^ (jKennicutt et all 

|200| ). 

iFraternali et al.l (|200ll ) have separated the so-called 
anomalous (extra-planar) gas i n NGC 2403 from the d isc gas 
with a technique described in iFraternah erahl (|2002| ). The 
resulting velocity field of the extra-planar gas (here shown 
in Fig. [5l upper panel) shows a clear sign of radial inflow in 
that the kinematical minor axis is rotated counter-clockwise 
and not orthogonal to the major axis. We have used the 
same procedure to derive velocity fields of the gas above 
the plane in our model cubes. The results are shown in the 
lower panels of Fig. [5l Clearly the inclusion of gas accretion 
changes the direction of the clouds' flow from outflow (bot- 
tom panel, pure fountain) to inflow (middle panel) . Note the 
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Figure 5. Velocity field fo r the extra-plan ar (halo) gas of 

NGC 2403 (upper panel) (see IFraternali et al.l I2OOII) compared 
with the two models of fountain + accretion (middle panel) and 
fountain only (bottom panel) . The velocity fields are rotated with 
the major axis of the galaxy along the horizontal axis. The tilt in 
the kinematic minor axis is the signature of infiow/outflow. 



tilt of the minor axis in the model is also quantitatively very 
similar to the tilt seen in the data. 

The above results occur because the orbits of fountain 
particles are significantly changed by their interactions with 
infalling gas; in Fig. |6] the top left and middle left panels 
show the orbit of a representative particle that is launched 
at ~ 60kms^^ from R ~ 5kpc in the pure fountain model 
(top) and the model with infall (middle; notice the different 
horizontal scales in the two panels) . In general the inclusion 
of accretion makes fountain clouds fall back to the disc at 
radii that are roughly the same as (or lower than) the start- 
ing radii, rat her than at larger radii as in cla ssical fountain 
models (e.g. ^Colli ns. Beniamin fc Rand|[2002l ). The bottom 
two panels of Fig. [S] show how in the presence of infall the 
particle increases in mass and loses specific angular momen- 
tum as a result of interactions. 

The same accretion model (with the same accretion pa- 
rameter) is able to produce both the amount of lag of the 
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Figure 6. Trajectories of a representative particle in tlie po- 
tential of NGC 2403. The upper panels show the case of a pure 
fountain, the middle panels show the effect of interaction with ac- 
creting material. In the bottom panel we show the modification 
of the mass and angular momentum of the fountain cloud as a 
consequence of this interaction. 



extra-planar gas in NGC 891 and the inflow pattern required 
by NGC 2403. Qualitatively it is clear that when fountain 
particles interact with the accreting material, they both have 
to share their angular momentum and acquire an inward mo- 
tion. However, it is remarkable that these effects are made 
quantitatively correct by a single value of a. 



5 RELATION TO OTHER WORK 

Our investigation of the coupling between clouds in the 
HI halo and a corona led to the conclusion that near the 
plane the corona would rapidly approach corotation with 
the disc. Galaxies with r apidly spinning , stead y-state coro- 
nae were investigated bv lBarnabe et aiTl|2006l ). They found 
that such coronae, taken in isolation, can reproduce the ro- 
tation curves shown in Fig. |31 raising the possibility that 
the observed HI clouds are embedded in and locally at rest 
with respect to a corona of this type. Because the coronae 
described by Barnabe et al. enjoy considerable rotational 
support, they are not very hot: in most of the volume above 
the disc, 2 x 10** K < T < 2 x 10^ K. At these temperatures 
radiative cooling is at peak efficiency, and unless their densi- 
ties are made extremely low, the cooling times ~ 1 kpc above 
the disc are much shorter than the local dynamical times. 
Similarly, if densities are not very low, the coronae would be 
readily detectable in UV emission lines of ions such as O vi. 

A second model that has recently been used to fit the 
rotation curves a bove the plane of NGC 891 is based on cos- 
mological inflow (|Kaufmann et aLll2006l ). In this model the 
gas is initially at the virial temperature and spinning. In a 



nearly spherical region of radius ~ 10 kpc the gas rapidly 
cools and collapses into the centre. Gas with low angular 
momentum moves in parallel to the spin axis to take its 
place, causing the temperature to rise as one approaches 
the axis at flxed z. Over several Gyr the remaining gas cools 
onto an annular region of the disc. At this stage the cool- 
ing gas lags the rotation of the disc in a manner similar to 
that observed in NGC 89 1. In fact, cooling from the initial 
conditions postulated by iKaufmann et al] (|2006l ) has cre- 
ated a rapidly spinning corona similar to that discussed by 
iBarnabe et~al] ||2006|} . 



Kaufmann et al.1 (|2006l ) argue that the cooling corona 



will spawn clouds of neutral H I that constitute the observed 
HI halo of NGC 891. Their simulations provide some sup- 
port for this conjecture, and they argue that the proposal i s 
supported by the analytic work of lMaller fc BuUockl (|2006l ). 
Since the simulations are subject to numerical artifacts as- 
sociated with their fln ite resolution in mass and distance 
(|Kaufmann et aLll2007t ). the analytic arguments for the ex- 
istence of cool clouds are important. 

We believe that these analytic arguments are flawed 
because they ignore the ability of a gravitational h eld to 
preven t a stratified atmos phere from suffering the iFieldl 
lll965h thermal instab ility (IMalagoh. Rosner fc Bodolll987l : 
iBalbus fc Sokerlll989l )F] Physically, when an atmosphere is 
confined by a gravitational field (which can be taken to in- 
clude the pseudo-field associated with spin), the gas set- 
tles on a dynamical time into a stratified configuration, in 
which the material of lowest specific entropy lies at the low- 
est potential. In this configuration an overdense region is 
a region in which the surfaces of constant specific entropy 
have been perturbed upwards. Consequently it is the crest 
of an internal gravity wave. On a timescale that is only 
slightly longer than the dynamical time, the overdensity 
sinks from lack of buoyancy, overshoots its point of equi- 
librium, and is soon an underdense region. Consequently 
the basis of the Field instability - the increase in cooling 
rate with increased density - cancels out to first order as 
initially overdense regions oscillate vertically and are on av- 
erage neither more nor less dense than undistur bed mate- 
rial at their mean position l|Balbus fc Sokerlll989l '). There is 
no reason to suppose that hot gas that is cooling onto the 
disc of a galaxy is more liable to the Field instability than 
is the gas of cooling flows, which X-ray spectra sho wed do 
not e volve according to the Field ins t ability (Peterso n et al.l 
I2002D . as iMalagoli. Rosner fc Bodd (Il987l ) predicted more 
than a decade earlier. 

A second problem with the contention that HI haloes 
represent material that is condensing out of the corona is the 
massiveness of haloes. In NGC 891, the HI mass detected 
above 1 kpc from t he plane is 1.2 x 10^ Afa, almost 30% of 
the total HI mass jOosterloo. Fraternali fc Sancisill2007f ). If 
this gas were falling towards the disk even at a relatively 
low speed of ~ lOOkms"^, taking a mean distance from the 
plane of about 4 kpc, it would imply an accretion rate of 
--^ 3OAf0yr^^. Not only does this rate substantially exceed 



1 In lMaller fc BuUockl r2006h the analysis relies on lBalbusI lll986t) . 
However in their equation (20) they misquote Balbus's equation 
(5) by misidentifying £(T, p) as the conventional cooling function 
A(T). 
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the current SFR, but it implies that we live in a time when 
the disc is growing much faster than it has in the past, which 
is a priori implausible and probably inconsistent with the 

measured col ours of the stellar d isc. 

Thirdly, iRand fc BeniaminI (120081 ) searched for HI 
around the only spiral galaxy that is known to have a large 
corona, the edge-on galaxy NGC 5746. Their finding that 
there is very little H I in the corona is hard to reconcile with 
the contention of iKaufmann et"al ] (HOOS) that thermal in- 
stability of coronae leads to the abundant formation of H I 
clouds. 

Finally, from an observational point of view there is 
strong evidence that a large fraction of the extra-planar gas 
is indeed produced by the galactic fountain. The distribu- 
tion of the extra-planar If I in NGC 891 and NGC 2403 is 
concentrated very close to the star-forming disk and almost 
absent above a nd below the ou ter disk. In other galaxies 
like NGC 6946 (|Boomsmal lioorl ) most of the high-velocity 
featu res are locat e d in the inner star forming disc. More- 
over, iHeald et all (|2007l ) found that in three galaxies with 
extra-planar ionised gas, the thickness of the ionised halo as 
well as the vertical rotational gradient both clearly correlate 
with the SFRs, pointing to a dominance of the fountain in 
the formation of these haloes. 



6 DISCUSSION AND CONCLUSIONS 

Standard cosmology teaches that the majority of the 
baryons lie somewhere in extragalactic space, probably in 
a hot diffuse medium. We cannot hope to have an adequate 
knowledge of the formation and evolution of galaxies until 
we have a better understanding of the connection between 
star-forming galaxies and extragalactic gas. In this paper 
we have shown that the existence of H I haloes around star- 
forming galaxies has profound implications for this vital con- 
nection. 

It has long been evident that star-forming galaxies 
must be accreting gas at significant rates. We have mod- 
ified the ballistic fountain model to include the effects of 
clouds sweeping up material that has low angular momen- 
tum about the disc's spin axis. The key parameter of the 
model is the rate a at which the mass of a cloud exponen- 
tiates as a result of accretion. We find that a single value, 
a « 1.5Gyr~^ enables us to resolve the problems encoun- 
tered in FB06 with both the edge-on galaxy NGC 891 and 
the inclined galaxy NGC 2403. Specifically, adding swept- 
up gas doubles the lag in rotation of the halo with respect 
to the disc because clouds have to share their angular mo- 
mentum with swept-up gas, and the any inward motion of 
accreted material causes the net motion of the fountain gas 
to be inwards rather than outwards. The aggregate infall 
rates implied for NGC 891 and NGC 2403 are ~ 2.9 and 
~ O.8M0yr~^, respectively. These rates are very close to 
the rates at which star formation is consuming gas in the two 
galaxies, but 1/5 — 1/10 of the rates at which star- formation 
cycles gas through the halo (FB06). This means that most 
of the neutral gas that we observe in the halo (~90% in 
NGC 891 and ~ 80 - 85% in NGC 2403) is gas that steUar 
feedback has pushed up from the disc and only a small frac- 
tion is accreted. Notably, the obtained accretion rates are 
not only consistent with the respective SFRs, but also with 



the expected cosmological accretion rates for these types of 
galaxies (e.g. iNeistein fc Dekelllioosi ). 

It is unclear whether the swept-up gas arises from cool- 
ing of a hot medium, or fr o m the infall of cold streams of gas 
jBirnboim fc Dekell l2003l : iBinnevI |2004 iKeres et all l2005l : 
ICattaneo et al.l 20061 ). but the long cooling time of coro- 
nal gas strongly suggests that it comes from cold streams. 
The accretion rates onto nearby spiral galaxies determined 
here are cosmologically significant: if the two galaxies had 
accreted at their current rates for a Hubble time, they would 
have added ~ 1/2—1/3 of their current masses. Hence our es- 
timates are consistent with past accretion rates being larger, 
but not much larger unless SN feedback is efficient in eject- 
ing material back to the IGM. Note that estimated accre- 
tion rates could be considered as lower limits because gas 
can reach the discs without encountering fountain clouds, 
either by slipping past these clouds if they have a small fill- 
ing factor, or by joining the disc at large radii, outside the 
star- forming region. 

The accretion rates estimated here exclude stars 
brought in by satellites, but since accreted stars cannot con- 
tribute to the thin disc, the bulk of the accretion onto these 
disc-dominated galaxies must be gaseous, consistent with 
the significant rates of gas accretion that we require. 

We have argued that feedback from star formation is 
essential for the formation of H I haloes - the latter cannot 
form directly from cooling coronal gas because (i) the latter 
is not thermally unstable, (ii) if H I haloes represented mate- 
rial that has condensed from the corona, the accretion rate 
of the thin disc would be implausibly large, and (iii) there 
is broad observational evidence of a link between galactic 
fountains and extra-planar gas features. 

In the likely event that coronal gas is not swept up by 
fountain clouds, the corona must be profoundly infiuenced 
by interaction with fountain gas. In particular, just above 
the disc interaction with halo clouds will ensure that the 
corona is nearly corotating. If the nearly corotating gas is 
confined near the disc, it must be quite cool (T ~ 10^ K) 
and strongly radiating in the UV. It seems more likely that 
the corona is hotter and therefore not bound to the vicinity 
of the disc: then it will flow outwards and upwards. 

We have to expect that the interface between HI in 
a galactic disc and coronal gas is at least as complex and 
dynamic as that between the solar photosphere and corona 
because similar physics is at work: global rotation, local vor- 
tices, a steep density gradient, upswelling of material, mag- 
netic loops and optically-thin radiative cooling in the tem- 
perature range 10* — lO'^ K. It is worth noting that above 
the Sun the transition from T = 3 x 10"^ K to T = 3 x lO'' K 
is complete within only ~ 30 km (|jordanlll977l ). 

Our assumption that swept-up gas has low angu- 
lar momentum about the disc's spin axis is both essen- 
tial and at first sight surprising. Indeed, discs must be 
formed by gaseous infall, and they are expected to form 
from inside out as a result of the mean angular mo- 
mentum of infalling material rising over time. Galaxies 
such as NGC 4550, which has two counter-rotating stel- 
lar discs l|Rubin. Graham fc Kennevlll992 ) and the numer- 
ous decoupled gas eous/stellar systems found with SAURON 
(jSarzi et al.l |2006| ) make it certain that whatever may be 
true of the time-evolution of the magnitude of the specific 
angular-momentum vector L of infalling gas, any given com- 
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ponent can diminish and even change sign. For our model 
it is essential only that at the radius R of accretion to a 
disc with circular speed Vc we have Lz < ^Rvc- The per- 
pendicular components of L could well make \L\ > Rvc- In 
the future, it will be possible to test these predictions with 
hydrodynamical cosmological simulations. Existing mod- 
els seem to agree q ualitatively with our requirements (e.g. 
IPorciani. Dekel fc Hoffmanll2002i l . 

Moreover, one interpretation of warps in HI discs re- 
lies precisely on L for infalling material being compara- 
ble in magnitud e to, but strongly misaligne d with, L for 
the outer disc llOstriker fc BinnevI [1989; Ji ang fc BinnevI 
1 19991 : IShen fc SellwoodI l2006h . This argument indicates 
that there could be a connection between disc warp- 
ing and rotational lag of the HI halo. Interestingly, 
both NGC 891 a nd NGC2403 do show warp ir ig of 
the outer discs (|Garci'a-Ruiz. Sancisi fc KuiikenI l2002l : 
iFraternali et al.ll2002l ). However, if the specific angular mo- 
mentum of infalling material is large enough, there could be 
negligible infall in the vicinity of the fountain gas. Then the 
galaxy would have a warped H I disc but an H I halo with 
only a small rotational lag. 

We have modelled the ambient gas as smooth, so foun- 
tain clouds receive a constant amount of mass and momen- 
tum per given time. In reality, an accretion flow is probably 
clumpy. However, as long as the gas mass of the accret- 
ing clouds is lower than the mass resolution in the data 
(about Mq), no significant difference is expected. For 
larger accreted masses, the effects of accretion will be more 
pronounced in some regions of the halo than in others, un- 
less large gas complexes break into smaller clumps, which 
can be modelled by smooth accretion. 

It would be interesting to extend this work to other 
galaxies, in particular low surface-brightness (LSB) galaxies 
with low st ar formation rates that also have HI haloes - 
UGC7321 l|Uson fc Matthew j|200j ) is an example of such 
a system. This would allow us to establish if accretion is a 
common process at the current epoch, and what role it plays 
in triggering star formation. We are also applying the model 
to the Milky Way; this work will clarify the relation between 
fountain gas and the classical phenomena of intermediate- 
and high-velocity clouds (Fraternali fc Binney, in prep.). 
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